Tilings of space and superhomogeneous point processes 
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Abstract 

We consider the construction of point processes from tilings, with equal volume tiles, of d- dimensional 
Euclidean space R d . We show that one can generate, with simple algorithms ascribing one or more 
points to each tile, point processes which are "superhomogeneous" (or "hyperuniform" ) , i.e., for 
which the structure factor S(k) vanishes when the wavenumber k tends to zero. The exponent of 
the leading small-fc behavior, S(k — » 0) oc fc 7 , depends in a simple manner on the nature of the 
correlation properties of the specific tiling and on the conservation of the mass moments of the tiles. 
Assigning one point to the center of mass of each tile gives the exponent 7 = 4 for any tiling in 
which the shapes and orientations of the tiles are short-range correlated. Smaller exponents, in the 
range 4 — d < 7 < 4 (and thus always superhomogeneous for d < 4), may be obtained in the case 
that the latter quantities have long-range correlations. Assigning more than one point to each tile 
in an appropriate way, we show that one can obtain arbitrarily higher exponents in both cases. We 
illustrate our results with explicit constructions using known deterministic tilings, as well as some 
simple stochastic tilings for which we can calculate S(k) exactly. Our results provide, we believe, 
the first explicit analytical construction of point processes with 7 > 4. Applications to condensed 
matter physics, and also to cosmology, are briefly discussed. 

PACS numbers: 02.50.-r, 61.43.-j, 98.80.-k 



I. INTRODUCTION 

"Superhomogeneous" [l] or "hyperuniform" Q point 
patterns in ci-dimensional Euclidean space M d are defined 
to be those in which infinite wavelength density fluctu- 
ations vanish. In other words, the structure factor (or 
power spectrum) S(k) of the number density field at wave 
vector k has the following behavior: 



lim S(k) = 0. 



(1) 



This defining characteristic of superhomogeneity (or hy- 
peruniformity) is tantamount to saying that the usual 
mean-square particle-number fluctuations increases less 
rapidly than R d for large R, where R denotes the lin- 
ear size of an observation window in R d [HQ. Indeed, 
the magnitude of such local density fluctuations have 
been suggested as a possible "order metric" to quan- 
tify the degree of order (disorder) of an arbitrary point 
pattern (2j. Any superhomogeneous point pattern can 
be seen as a typical configuration of a particular type 
of "critical" point in that the direct correlation function 
(defined through the Ornstein-Zcrnikc relation) is long- 
ranged while the pair correlation function is short-ranged 
0. Such remarkable behavior is diametric to that seen 
in usual thermal critical points in which the inverse is 
true, i.e., the pair correlation is long-ranged and the di- 
rect correlation function is short-ranged. 



Although it is clear that any periodic point pattern 
is superhomogeneous, it is less obvious that statistically 
translationally and even rotationally invariant random 
point patterns in M. d can have this property. We now 
know of a variety of intriguing translationally and ro- 
tationally invariant random point patterns that are su- 
perhomogeneous, including the ground state of liquid 
4 He [H, lilfEj], maximally random jammed hard-sphere 
packings U certain one-component plasmas 0, H, 0]j 
the matter distribution in the Universe [l|, [8|, and cer- 
tain aperiodic tilings 0, H, [l(| . An interesting applica- 
tion of superhomogeneous point patterns in cosmology is 
in the preparation of initial conditions for gravitational 
A-body simulations [1, [^, [l(| EH- Superhomogeneous 
distributions also appear in cosmology in the context of 
the determination of bounds, first derived by Zeldovich 
12], on the mass fluctuations at large scales generated 
by causal mechanisms (i.e. with physics respecting the 
causal constraints of cosmological models). Indeed, we 
note that in this context a simplified form of the analysis 
we develop here of the small-A; behavior of the structure 
factor is often used (see e.g. Ref. [131]). 

It is desirable to develop both theoretical and compu- 
tational methods to generate a wide class of superhomo- 
geneous random point patterns. Recently, a collective 
coordinate approach [3, [l5| has been employed to nu- 
merically generate translationally invariant superhomo- 
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geneous point processes. This procedure enables one to 
produce point patterns that completely suppress density 
fluctuations of modes for a positive range of wavenum- 
bers around the origin. In Ref. [l6| an algorithm for gen- 
erating discrete processes in one dimension with super- 
homogeneous mass fluctuations has been given (see also 
Ref. [IZl)- An analytical methodology to relate super ho- 
mogeneous point processes and Voronoi tilings of space 
has recently been proposed and studied in Ref. (l8| . 

In this paper, we study the construction of superho- 
mogeneous point patterns starting from generic tilings 
of Euclidean space R d with equal volume tiles. We 
show how to explicitly generate such point processes in 
which the structure factor for small wavenumbers has the 
power-law form S(k) ~ fc 7 for positive 7, where fc = |k| 
is the wavenumber. The constructions illustrate the very 
specific properties of these superhomogeneous point pat- 
terns in which the exponent 7, characterizing the long 
wavelength fluctuation in fc space, is related to the de- 
tailed arrangement of the points on small scales. Our 
study also shows how the exponents of the small-fc be- 
havior of the structure factor for these point processes 
encode properties of the tilings, and could thus possibly 
be used as a method for classifying them. In a related 
article by two of us [l9[ the two-point correlation proper- 
ties of point processes generated by replacing each parti- 
cle, in a point process with known two point properties, 
by a "cloud" of particles are derived 1 . One of us [15| nu- 
merically generated disordered point distributions within 
a cubical box under periodic boundary conditions with 
7 > 4. However, to our knowledge, prior to this paper 
and [l9l ]. explicit analytical constructions of point pro- 
cesses with 7 > 4 have not been given previously in the 
literature. 

It is instructive to recall qualitatively why tilings are 
a natural starting point for the construction of superho- 
mogeneous point processes. A tiling or tessellation is a 
partition of Euclidean space M. d into closed regions whose 
interiors are disjoint regions [2fl| ■ Let us suppose we have 
a tiling of space by tiles which are (i) of equal volume 
||T||, and (ii) bounded, with maximal length A in any 
direction. Let us now place one point in each tile and 
consider the number fluctuations in the point process so 
generated. If N(R) is the number of points in a sphere 
of radius R, and of volume V(R), it is simple to see that 

ViR - A) <N(R)< ViR+ u A) . (2) 



1 Results in this case are derived in |19H under the assumption that 
the stochastic process describing the generation of the "clouds" 
and the initial point process are independent. In the algorithm 
discussed here this is not the case, as the points are ascribed to 
each tile in a way which depends, in general, on the tile. For the 
particular case of a Bravais lattice tiling, however, both calcu- 
lations are valid because of the equivalence of all tiles/points in 
such a lattice. Indeed, in this case, the different general formulae 
derived in the present article and [l9ll give the same result. 



The lower bound is the minimal number of tiles which 
can overlap the sphere of radius R — A [and all such tiles 
must contribute a point to N(R)}, the upper bound is the 
maximal number of tiles which are fully enclosed in the 
sphere of radius R + A [and only such tiles can contribute 
to N(R)]. For R — ► 00 we have therefore 



\AN(R)\ < cR d ~ 



(3) 



where c is a constant, and AN(R) = N(R) - N(R) with 
N(R) = V(R)/\\T\\. Averaging over configurations (or 
randomly placed centers for the spheres) one anticipates 
that the slowest possible scaling of number fluctuations is 



(AN 2 {R)) cx R 



d-l 



(4) 



where (...) denotes the ensemble (or volume) average. 
This behavior of the variance, proportional to the sur- 
face, is a characteristic of superhomogeneous point pro- 
cesses. If there is appropriate long-range correlation in 
the tiling at arbitrarily large scale, the fluctuations could 
however, in principle, add coherently to give the more 
rapid growth up to 



(AAf 2 (i?)) cxi? 2 ^ 1 ', 



(5) 



which corresponds to the limit equality in Eq. . While 
in d — 1 this still corresponds to surface fluctuations 2 , 
for any d > 2 it implies only the limiting small- fc behav- 
ior S(k) oc fc 7 with 7 > — d + 2, which means that the 
point processes are not necessarily superhomogeneous for 
d > 2. We will recover this result below, with the only 
difference that the bound is found to be 7 > — d + 4, 
which implies that even long-range correlations between 
tiles give superhomogeneous processes for d < 4. The 
difference between this result and our naive estimate is 
simply due to the fact that below we constrain the par- 
ticles to lie at the center of mass, rather than placing 
them randomly. In fact we will show here that by assign- 
ing more than one point in an appropriately constrained 
manner to each tile, we can increase these bounds on the 
exponents without limit, and realize superhomogeneous 
processes with an arbitrary positive exponent, in any di- 
mension. 



II. POINT DISTRIBUTIONS FROM TILINGS: 
ONE POINT PER TILE 

We consider in this section point processes generated 
by ascribing one point to each tile. We first give a general 
analysis of the small-fc properties of the structure factor 
of the density fluctuations, and derive how the leading 
behavior is determined by the properties of the tiling. We 
then describe some specific explicit constructions which 
illustrate the result. 



2 Note that the case of d = 1 is rather trivial given our assump- 
tions: the only equal volume tiling is the lattice. 
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A. Density fluctuations and long-wavelength limit 

We start from a generic (regular or irregular) tiling 3 of 
rf-dimcnsional Euclidean space M. d into equal volume tiles, 
which we denote T^. We consider the point distribution 
generated by ascribing one point to each tile, and placing 
it at position x 1 , which coincides with the center of mass 
of the tile T,-, i.e., 



1 

iirii 



d d a 



(6) 



where ||T|| is the volume of the tiles. The density fluctu- 
ation field is thus 



n 



(7) 



where no is the mean number density in the infinite- 
volume limit 4 , and <5^(x) is the Dirac delta function 
in d dimensions. The structure factor (SF) is defined as 



3(k)= lim feZ^ = l + noMk) 

V— >oo UqV 



where V is the system volume, 
in(k; V) = 

and 



d d xe" lk x Sn(x) , 



h(k) = / d d 



-ik-x 



Mr) 



(8) 



(9) 



(10) 



is the infinite-space Fourier transform of the total pair 
correlation function h(r), which vanishes for disordered 
systems when the distance r = |r| tends to infinity 
For our point process it follows directly that 



<fn(k; V) = e "' k ' Xl l 1 ~ ^( k ) 



(11) 



where the sum runs over the points enclosed in the vol- 
ume V, and Wj(k) is the normalized characteristic func- 
tion of the tile Ti, given by 



Wi(k) 



1 

iirii 



d d xe~ 



(12) 



Ti(0) 



where Ti(0) denotes that the center of mass of the tile 
has been taken as the origin of axes. If we assume that 



3 A regular tiling is periodic in space. An irregular tiling is ape- 
riodic in space, including quasiperiodic as well as disordered 
tilings. A congruent tiling consists of identical tiles. 

4 Since all particles have the same mass no distinction need be 
made between the mass and number density fluctuations. For 
the case of a single particle per tile no = 1/||T||. 



Wi is an analytic function at k = 0, we can expand it in 
Taylor series, to obtain 



^— ' ml 

m=2 



M ai ... am (i) (13) 



where 



M ai ... c 



>(*) - || T |, 



llj 00 fy r, m m m 30 Q 



(14) 



Ti(0) 



is a (fully symmetric) tensor of rank m corresponding to 
the m-th moment of the mass distribution of the tile T, 
(normalized by the volume/total mass) and ctj = 1, . . . , d 
are indices for the Cartesian components. Note that we 
have used Eq. ([6]), which makes the linear term in the ex- 
pansion (|13p (corresponding to the dipole moment) van- 
ish. The assumption of analyticity corresponds to the 
requirement that all these moments are finite. This is 
true in particular if the tiles are of finite extent. We 
will discuss briefly in our conclusions the possibility of 
relaxing this assumption. 

Using these expressions in the definition of S(k) we 
now obtain 



S(k) 



(15) 



(—i) m (i) n 

T~i k ai ...k an k/3 1 ...ki3 m l ai ... an [3 1 ... l 3 m ('k) 

n=2m=2 

where 

T ai ... an01 ... 0m (k) = ^lim iEE e - ik - (X< ^' } 

i 3 

xM ai ... an (i)M l3l ... l 3 m (j) (16) 

where the sums run over the TV particles contained in 
the volume V. It is straightforward to verify that the 
coefficient of the leading term in k (at order k A ) is non- 
negative, and that the coefficients of all powers of k are 
real. Indeed the SF 5(k) is by definition a real non- 
negative quantity, and Eq. (|15p is just the specific form 
of its Taylor expansion around k = for the particular 
class of distributions we are considering. 

It is convenient to rewrite the latter expression as 

r Ql ...a„/ 3l .../3 m (k) = lim J- [ d d xd d ye- ik <*-^ 

M W^ft... fc (y)(i7) 

where 

M Oil . . .Oi n 

(x)=V^)(x-x l )¥ a ,.. a „( ! ). (18) 



The distribution M. (x) can be viewed as a weighted par- 
ticle density. The weight associated with each particle is 
the appropriate component of the mass moment of the 
tile to which the particle belongs. 
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Up to now we have considered implicitly a single parti- 
cle placed deterministically in each tile. We now consider 
averaging over an appropriately defined ensemble of such 
tilings 5 . If the tiling is statistically translationally invari- 
ant, we have that 

(M ai ... an (^)M l3l ... 0m (y)) = 3 "< m (x - y) , (19) 

where (...) denotes the ensemble average. We have 
adopted here for the correlation function g"' m (x) the ten- 
sorial notation in which the indices are left implicit. In 
this notation we can write our result for the SF S(k) as 

oo oo t_j\mfj\n 

s ( k ) = -EE^r k "'Hk)'k™ (20) 

U n=2 m=2 

where <?™' m (k) is the Fourier transform of g" ,m (x) defined 
as 



.9" 



l (k) = J Ae- ,k ' x f- m (x) 



(21) 



and k" denotes a tensor of order n, given by the tensor 
product of n vectors k, i.e., 

k ai ...a n — [k®k-- ®k]ai...a„ — k a2 ■■k ctn . (22) 

The symbol • in Eq. (f2T)]) denotes the contraction of the 
corresponding tensor indices. If the ensemble is also sta- 
tistically isotropic the product k™ • <7™ ,m (k) -k m , and thus 
S(k), is a function of k — |k| only. 

The behavior of <S*(k) at small k is thus manifestly de- 
termined by that of g n,m (k) in this limit. These quanti- 
ties are in fact the (tensor) SFs associated to the discrete 
stochastic field defined by Eq. (| 18[> . the two point cor- 
relation function of which is g n,m {x). They thus encode 
information about the tiling, and more specifically about 
the correlation properties of the second and higher mo- 
ments of the tiles. We restrict ourselves to the case that 
all these moments are finite, and strictly bounded (which 
also ensures, as noted above, the validity of the expansion 
of the characteristic function Wi (k) we have performed) . 
As noted above, each component of the tensor stochastic 
field M n (x) defined in Eq. (fT5|) above is then a discrete 
stochastic process in which the points are located at the 
same positions as in the point process we are studying, 
but have "masses" given by the corresponding component 
of the tensor M n {i) (or rather "charges" as they are not 
strictly positive) which are bounded (above and below). 
Just as for a generic stochastic point process such a dis- 
crete (or indeed continuous) process can be classified into 
three categories according to the small-fc behavior of the 



1. g n ' m (k = 0) = const. > 0: this means that the cor- 
relation functions g™' m (x) of the higher order mo- 
ments are integrable at large x, and the correspond- 
ing integral is equal to a positive constant, i.e., the 
higher moments of the tiles have short-range cor- 
relations dominated by the positive contributions. 
For the generated point process we have then the 
leading behavior S(k — > 0) oc k 4 . 

2. <7™' m (k = 0) = 0: the integral of the correlation 
functions g n ' m (x) converges to zero, i.e., the shapes 
and orientations of the tiles have themselves super- 
homogeneous properties (i.e. in which the positive 
and negative correlations balance exactly in the in- 
tegral) . In this case we will obtain a leading behav- 
ior S^k — > 0) oc fc 7 with 7 > 4 [and with a value 
depending on the leading behavior of the g n ' m (k) 
at small- fc]. 

3. g n ' m (k = 0) = oo, with <?™> m (k -> 0) oc k a and 
— d < a < 0. In this case, in which the correla- 
tion functions g n ' m (x) are non-integrable, i.e., the 
higher moments of the tiles have themselves long- 
range correlations, we can obtain a leading behav- 
ior for our point process 5(k — > 0) oc fc 7 with 
4 - d < 7 < 4. 

Several remarks on this result are important. Firstly, 
for simplicity we have assumed above that the g n ' m (k) 
are in the same class for all m and n. This is, of course, a 
priori, not necessarily the case. In the more general case 
that the different <7™' m (k) are in different classes, the de- 
termination using Eq. (|20[) of the exponent 7 of the lead- 
ing small-fc behavior of S(k) is nevertheless straightfor- 
ward. Furthermore, it is simple to verify that the bounds 
we have given on this exponent remain valid. 

Secondly, we have assumed that the g n ' m (k) obey the 
condition 



lim fcV 



'(k) = 0. 



(23) 



This assumption corresponds to the requirement that the 
discrete processes defined by Eq. (IT51) have well-defined 
mean values, i.e, the normalized fluctuations (e.g. inte- 
grated in a sphere) of the moments of the tiles converge 
to zero in the infinite- volume limit. While this seems a 
very weak assumption, it is not a priori true of all tilings. 



B. Explicit constructions 

We now give various explicit constructions to illustrate 
the above results. 



5 For a deterministic tiling, e.g. the cells of a regular lattice, or the 
pinwheel tiling in two dimensions discussed below, this average 
can be defined by the set of configurations generated by applying 
an arbitrary rigid translation to a given configuration. 



1. Regular lattice tilings 

Consider first the tiling given by the Voronoi cells of 
any Bravais lattice. In a Bravais lattice, the space 3t d can 
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be geometrically divided into identical regions F called 
fundamental cells, each of which contains just one point 
of the lattice. A Voronoi cell associated with a point at 
r in a point distribution is defined to be the region of 
space nearer to the point at r than to any other point. 
Since all Voronoi cells or tiles of a Bravais lattice have the 
same shape and orientation, the mass moments M n (i) in 
Eq. ([H]), calculated with respect to the center of mass, 
are identical for all tiles, i.e., M n (i) does not depend on 
i. The discrete processes specified by Eq. (fig]) are then 
simply, up to a constant, equal to the density field of 
the lattice, and thus it follows that all the correlation 
functions g n,m (x) are proportional to the two point cor- 
relation function of the original lattice. Thus <?™' m (k), 
and also S(k) = 0, is zero in some finite region around 
k = 6 . This result is in fact evident: the point process 
generated by placing points at the center of mass of every 
cell is of course simply the lattice itself. In this case, of 
course, neither the tiling nor the point process are sta- 
tistically isotropic (where the statistical average is taken 
over lattices rigidly translated within the elementary lat- 
tice cell) 

The same result can evidently be generalized to any 
tiling of equal volume cells constructed from a periodic 
point pattern. 



2. Congruent rotationally invariant tilings 

We consider next deterministic congruent tilings with 
the additional property of rotational invariance, i.e., in 
which all orientations of the identical tiles are equiproba- 
ble. Known examples are the pinwheel tiling |21j in d = 2 
(see Fig. 1) and the quaquaversal tiling [22j in d = 3. In 
these cases each tile Ti can be characterized solely by its 
center of mass Xj and by a matrix R(i) £ 0(d), the lat- 
ter giving the orientation with respect to some arbitrary 
chosen orientation. 

We can then write, in tensorial notation, 

M n (x) = ^<5 (<i) (x - Xi)iT(i) • M n (24) 

i 

where R n (i) is the tensor product, giving a tensor of rank 
2n of which n indices are contracted with the correspond- 
ing moment Ai n of a tile with the reference orientation, 
i.e., 

[R(i) n ■ M n ] ai ... an = R ai01 {i)...R anPn {i)M n Pl ... 0n (25) 

where, as everywhere above, the sums over the indices 
which appear twice are implicit. The correlation func- 
tions g n ' m (x) are thus direct measures of the correlation 




FIG. 1: Portion of a pinwheel tiling. The prototile of the 
pinwheel tiling is a right triangle with sides of length one, 
two, and vo. The tiling is produced by performing certain 
"decomposition" and "inflation" operations on the prototile. 
In the first step, the prototile is subdivided into five copies of 
itself and then these new triangles are expanded to the size of 
the original triangle. These decomposition and inflation oper- 
ations are repeated ad infinitum until the triangles completely 
cover the plane. 



of the orientations of tiles with center of mass separated 
by x. 

The case of the quaquaversal tiling has been studied 
numerically in Ref. lOj , and an approximate small-fc be- 
havior S(k) oc k 4 found. In Ref. [23| it is noted, however, 
that the numerical results agree better with 5(k) oc k 1 
and 7 rs 3.4. While the former behavior would corre- 
spond, as discussed above, to a short-range correlation 
of the orientation of the tiles, the latter would instead 
correspond to a weak long-range correlation (with a cor- 
relation function characterizing the orientations decay- 
ing with distance r as ~ r~ 2A ). Further, superimposed 
on this power-law behavior there are residual peaks at 
certain wavenumbers, with a spacing which appears to 
be consistent |10f with the hierarchical nature of the 
tiling [24|. As we have discussed the small-fc behavior 
of the constructed point process thus probes the correla- 
tion properties of the underlying tiling. 



3. Random Binary Rectangular (RBR) tiling 



More precisely 5(k) = at all k different from non-zero recipro- 
cal lattice vectors. 



It is instructive to illustrate our result with a non- 
trivial example which, albeit not statistically isotropic, 
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allows us to calculate exactly the SF of a point process 
constructed by the algorithm we have described. The 
example we now give is of a stochastic congruent tiling. 
For simplicity we work in d = 2, but a generalization to 
any d is straightforward. 

We generate the tiling as follows. We start from a 
regular tiling of the plane with congruent squares. We 
then divide each square tile in half, defining two identical 
rectangular sub-tiles, as shown in Fig. [2 The choice 



1 



m=l 




i 


• 


• 


i 



underlying square lattice is I. Consequently the volume of 
the elementary rectangular tile is l 2 /2 and therefore the 
average number density of the point process is no = 2/1 2 . 
It is simple to show that the Fourier transform (FT) of 
[n(x) -n ] is 



Sn(k; V) = 2 ^ e" 



ik-R 



R 



cos k~l 



,1 



- h kyl — 



OR 



six\(k x l/2) sm(k y l/2) 



kx^ /2 kyl i ^2 



(27) 



To calculate the SF averaged over the ensemble of pos- 
sible configurations of the binary tiles we assume that 
(<tr) = (i.e. both orientations of the binary tiles 
are equiprobable), and write x(R) = (or <7r +r) using 
the lattice statistical translational invariance. Moreover, 
since ctr = ±1, we have that 



cos k T l 



,1 + tr 



8 



1 + CTR 



4 



1 °R 



1 - CR 



(28) 



G=-l 



FIG. 2: Elementary binary rectangular tiling, with a unit 
mass particle in the center of mass of each rectangle, and its 
description in terms of an Ising-like spin variable. 



of the orientation of each tile is given by a stochastic 
process, which can be cast as the value of a simple up- 
down spin variable. The density of the point process 
generated using the algorithm analyzed in the previous 
section, in which a point is placed at the center of mass 
of each tile, can then be written 



It is then straightforward to obtain the following exact 
expression for the SF: 



S(k) 



where 



[cos(M/4) - cos(fc^/4)] 2 



Sml (k) 



R 



x(R) 

(29) 



S A fL(k)=A(k)^e- lk - R 

R 



7T 2 J2 A(H)S(k-H) 

H#0 



(30) 

is a "modulated lattice" SF which is different from zero 
only at the non-zero reciprocal lattice vectors H, and 



A(k) = \ 



cos 



4 



cos 



|Sin(fc x //2) sin(fcj,Z/2) 



(31) 



At small-A; only the first term of Eq. (f29[) contributes, 
and expanding the factor outside the sum we thus obtain 
a leading small- k behavior 



>(*) = EE* 



R n=l 



xS 



y- Ry 



2 

1 - O-R 



(-1)' 



(26) 



where R = (R x , R v ) are the lattice sites of the underlying 
square lattice placed at the center of each square cell, and 
(Jr = ±1 is the spin variable specifying the orientation 
of the two elementary rectangular tiles at the lattice site 
R as in Fig. [5] We assume that the lattice spacing of the 



S(k) 



where 



y(k) = hm — 

R 



,— ik-R 



x(R). 



(32) 



(33) 



This exact result is of course a special case of the general 
analysis given above, in which the moments characteriz- 
ing the tiles are particularly simple as there are only two 
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orientations. All the two-point properties of the point 
process are then contained in the single correlation func- 
tion x(R) of these orientations for pairs of tiles with cen- 
ters separated by R. As x(k) is a power spectrum of a 
stochastic process with a well defined mean, at small fc 
we have x(k) ~ k b with b > -d = -2. If x(O) = c > 0, 
i.e., in the case in which the orientations of the tiles are 
short-range correlated, then S(k) — (3(fc 4 ) at small k. 
If instead \(k — > oo) = oo there are long-range positive 
correlations in the orientations of the tiles which induce 
a slower decay of the fluctuations in the associated point 
process at large scales. Finally if x(0) = the stochastic 
spin process is itself supcrhomogeneous, with a balance 
between positive and negative correlations creating a sort 
of "stochastic order" in the spin configuration. In this 
case S(k) vanishes faster than fc 4 at small k. In Appendix 
|X]we describe explicitly an algorithm for generating such 
spin configurations. 



III. 



POINT DISTRIBUTIONS FROM TILINGS: 
n > 1 POINTS PER TILE 



The algorithm described in Sect. Ill B 31 can in fact be 
thought of in a different way to that in which we have 
presented it: one can consider it instead as a direct as- 
signment of two points to the tiles of the original square 
lattice, without any construction of an intermediate SBR 
tiling. The pair of particles then have two possible ori- 
entations, which are chosen stochastically. Note that in 
each case the center of mass of the particles is located 
at the center of the square cell. Our results above shows 
that if this stochastic process is short-ranged correlated 
we obtain a small-fc behavior S(k) oc fc 4 , while with a 
single point at the center of mass of the lattice cell we 
recovered (evidently) the lattice. We now consider quite 
generally what small-fc behavior of S(k) of a point pro- 
cess we can obtain by ascribing more than one point per 
tile in a generic tiling with equal volume tiles. 



A. Density fluctuations and long-wavelength limit 

We ascribe p points to each tile, denoting their posi- 
tions by 



x 4 ^ = x l 



u ie (34) 

with £ — l...p and x 1 is, as above, the center of mass 
of the tile 7$ (so that u 1 ^ is the position relative to the 
center of mass). We assume further that the center of 
mass of the points coincides with that of the tile, i.e., 



0. 



(35) 



i=i 



Following the same steps as in Sect. HH we arrive at 



Sn(k; V) = e 



1 p 

-y, < k " 



(36) 



where Wj(k) is precisely the same normalized character- 
istic function of the tile as defined in Eq. (Tj"2"|) [and 1/p 
is the fraction of the mass of the tile ascribed to each 
particle] . Expanding in Taylor series we obtain 



P 



=i 

where now 



lyy w,£ -wi(k) = V ^-k ai ...k am M ai ... am (i) 

(37) 



M ai ... am (i) = - 



1^ 



P 



1=1 



is the totally symmetric rank m tensor corresponding to 
the difference of the m-th moments of the mass distri- 
bution of the points associated with the tile Tj and that 
of the tile itself. As in the derivation with one point, 
we have assumed the analyticity of the quantity we ex- 
panded. This means that we require that all the moments 
in Eq. I|38p are finite, which is true in particular if all the 
tiles are of finite extent and the lengths of the vectors 
u 1 ' are bounded. 

All the expressions, and notably the result for S(k) in 
Eq. (|2T))) . derived in the case with one point per tile, are 
thus valid. The only difference is that the tensors M n (i) 
[and correspondingly A4 n (x)] are now given by Eq. (|38[) . 
The small-fc properties thus depend, assuming statisti- 
cal translational invariance, on those of the FT of the 
correlation functions of the differences of the moments 
of the discrete mass distribution of the points ascribed 
to each tile and that of the continuous mass distribution 
represented by the tile itself. 

The most important implication of this result is the 
following: the coefficients in the small-fc expansion of 
Eq. ([57)) are now proportional to a difference of two quan- 
tities in Eq. (|38[) . Any given coefficient will vanish iden- 
tically if our assignment of the p points satisfies the con- 
straints 



^S u of-«^ 6 = liiir / d d xx ai ...x an 
Pf^i \\T\\J Tl (o) 



(39) 



in each tile, i.e., if the (tensorial) moments of the mass 
distribution of the points are equal to those of the tile in 
which they are placed. With a sufficiently large number of 
points per tile one can evidently make any desired finite 
number of terms vanish in the expansion of S'(k), i.e., 
construct a stochastic point mass distribution for each 
tile which has all moments up to a certain order equal to 
those of the continuous mass distribution of the tile. 

This procedure allows one to obtain an arbitrarily large 
exponent 7 in the small-fc behavior of the SF of the point 
process. 



B. Explicit constructions 

We again illustrate these results with some explicit ex- 
amples. 
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1. Regular lattice tiling 

In a Bravais lattice, as we have discussed above, the 
tensorial moments of all tiles are equal so that the second 
term in Eq. (|38|) does not contribute to S(k) in a finite 
region around k = 0. This can most easily be seen by us- 
ing Eq. (|38|) directly in the expressions Eqs. (fT5|) and (fT6|) 
for S'(k): the i-independent second term in each compo- 
nent of the tensor M n (i), when summed over i, gives a 
delta function proportional to the SF of the lattice. If we 
have more than one point per tile (i.e., p > 2) we may 
have, however, an i-dependent contribution from the first 
term in Eq. (|3"5)l . i.e., from the moments of the mass dis- 
tribution constituted by the points assigned to the single 
tile, which are not constrained (beyond the dipole mo- 
ment). If we allow the distribution of these points to 
vary stochastically from cell to cell, we will generically 
have a non-zero contribution for the SF (ensemble aver- 
aged over the stochastic process) at all k, i.e., we will 
have a continuous SF 7 . 

As a simple example let us consider first the RBR algo- 
rithm analyzed above, cast as the ascription of two points 
to each cell of a simple cubic lattice, but now allowing the 
pair of particles ascribed to each cubic cell have a ran- 
dom orientation, i.e., we take two points in each lattice 
cell with coordinates 

xi (R) = R + u(R) , x 2 (R) = R - u(R) (40) 

where R is the generic lattice site (which is the center 
of mass of the corresponding tile), and the vectors u(R) 
are generated by a stochastic process. We assume fur- 
ther that the vectors u(R) in different lattice cells are 
uncorrelated. The ensemble is thus fully specified by the 
one point probability distribution function p(u). We call 
this stochastic point process the split shuffled lattice 8 , 
as it is a generalization of the shuffled lattice discussed 
in Ref. [l| (see also Refs. [ill, HH), m which one point is 
randomly displaced off a perfect lattice. 

The leading small- fc behavior of S(k) in this case may 
be found easily by taking the ensemble average of the 
leading term in Eq. (TT5)) : 

S(k) = ifc a fc /3 fc 7 fc 5 (X Q/ 3 7< 5(k)) (41) 



If, on the other hand, the points are placed in the same way with 
respect to the center of mass of each tile, 5(k) is again zero in 
the same region. The remaining non-zero piece has a modulated 
delta-function structure which can be easily calculated. Indeed 
the point distribution so generated is in this case again a periodic 
lattice, i.e., the initial Bravais lattice with basis of p points per 
cell. 

In the context of the context of causality bounds on fluctuations 
in cosmology, this construction has been studied in Ref. |25|. 



where, using Eq. (j38j) in Eq. (Ti"6)) , we have 

(2a/s 7 *(k)> = lim ^ E E e-* k - (R - R,) 

R R' 

x (u a (R)up(~R.)u 7 (R')u6(R.')) . (42) 

Since the vectors u(R) are, by assumption, uncorrelated 
at different sites, we have 

(M ct (R)u /3 (R)u 7 (R')u ( 5(R')) = (u Q u / 3)(u 7 -u 5 ) . (43) 

for R 7^ R'. Using the fact, again, that the sum 
~k ErEr' e - lk -( R - R ) i s proportional to the SF of the 
original lattice, which is zero around k = 0, we can then 
write the leading small- k behavior as 

S(k) = l[((k.u)V((k.u) 2 ) 2 ] . (44) 

If we choose a probability distribution which is isotropic 
in u, i.e., p(u) = p{u) (with u = |u|), we then have 

S(k) = ^ [C(d)(u 4 ) - (uY] . (45) 

where C(d) > 1 is a constant [(7(1) = 1, C(3) =9/5]. 
We thus obtain the same exponent 7 as in the RBR tiling 
model, but now the fluctuations at small k are isotropic 
at leading order because of the isotropy in the attribution 
of the displacement vectors. Note that we have assumed 
here, as we do throughout this paper, that the moments 
of the mass distribution in each tile are finite, which re- 
quires here manifestly in Eq. (j45"]) the finitcncss of at least 
the first four moments of p{u). It is simple, however, to 
generalize this kind of model to the case when moments 
of order lower than the fourth diverge, just as done for 
the shuffled lattice in Ref. [26| (see also [Ty|). In this case 
one can obtain any small-fc behavior (0 < 7 < 4). 

This kind of algorithm can easily be generalized, con- 
serving a sufficiently large number of moments in such a 
way as to obtain higher powers of the small-/c behavior, in 
principle producing any desired leading behavior. Let us 
suppose that we generate mass moments M™ (R) at each 
lattice site R defined by Eq. (|3"5|) , with an uncorrelated 
stochastic process, i.e., so that 

(M"(R)M m (R)) - (M n (R)) (M m (R')) (46) 

for each tensor component of the tensor product, for R ^ 
R', (and any n, m). As in the example given above 
starting from Eqs. (fT5|) and (Tl"rj)) . it is straightforward to 
show that S'(k) may then be written, in a finite region 
around k = 0, as 

5(k) = E E tlL n L ( 47 ) 

n—2 rn — 2 

xk" • [(M n M m ) - (M n )(M rn )} ■ k m 

where all the tensors M n are evaluated at the same ar- 
bitrary lattice cell, i.e., 

[(M n M m )] ai ... an ^... 0m = (M ai ... an (R)M£ „ (R)) 

(48) 
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etc.. 

Thus, for example, we can obtain 7 = 6 with an algo- 
rithm of this type which allows the third moment M 3 of 
the mass distribution to vary from cell to cell while keep- 
ing the second (quadrupole) moment M 2 fixed. This can 
be done in a rotationally invariant manner by using a con- 
figuration of points with a quadrupole moment (relative 
to its center of mass) proportional to the identity matrix, 
e.g., points placed at the corners of a regular tetrahedron 
in d dimensions. Placing such a configuration at each lat- 
tice site, but rotated by a random rotation in SO(d), the 
variance in the second moment of the cell is thus zero. 
However it is easy to verify that the random rotation 
engenders a variation in the third moment, so that one 
obtains the leading small-fc behavior 5(k) cx fc 6 . The co- 
efficient of the fc 6 term can most easily be made zero by 
taking instead a configuration whose quadrupole moment 
is again diagonal, but whose third mass moment M 3 is 
zero. Since the latter is in fact zero for any configuration 
of points which is invariant under inversion symmetry, a 
possible choice is to add to each tetrahedron a reflection 
of itself in its center of mass. Alternatively, and using 
fewer points, one can use a randomly rotated configura- 
tion consisting of d couples of points with equal separa- 
tion placed orthogonally with common center of mass. 
In this case we obtain a leading order small-fc behavior 
5(k) cx fc 8 . Further discussion of this kind of point pro- 
cess generation can be found in Ref . . In this context 
they arise as a special case of "cloud processes" , in which 
points in a generic initial point process, of which the cor- 
relation properties are assumed known, are replaced by 
a "cloud" of point particles. 



2. Congruent rotationally invariant tilings 

We have seen that superhomogeneous point processes 
with a small-fc behavior 5(k) cx fc 7 and 7 > 4 can be 
generated starting from a regular lattice tiling, by using 
a stochastic process to determine the positions of an ap- 
propriately constrained set of points in each cell. While 
the point process so generated is not statistically transla- 
tionally and rotationally invariant, we saw that a leading 
behavior proportional to k was obtained if the stochastic 
process assigning the points had itself no preferred direc- 
tion. Thus at large scales the system approximates very 
well statistical translational and rotational invariance. 

Starting from a translation and rotation invariant 
tiling we can obtain a statistically translation and ro- 
tation invariant point process in an analogous way. It 
suffices in this case, however, to assign the points deter- 
ministically to each tile, in the same way relative to each 
tile (e.g. at the center of mass of the tile). We do not 
need now the additional stochasticity provided either by 
taking more points per tile or random positions inside the 
tile. The reason is that the moments of the tiles given 
by Eq. (TTJJ already depend on the tile itself because of 
the orientation. Consequently the quantities g" ,m (k) are 



non-zero around k = just as in the case when we had 
one point per tile. The difference is, as we have discussed 
at length above, that we can make certain terms zero so 
that the leading term appears at higher order in k. 

To see how this can be done in a little more detail, let 
us consider, to be specific, a pinwhcel tiling (d = 2) or 
quaquaversal tiling (d = 3) as in Sect. Ill B 21 above. It 
is convenient to study the SF 5(k) as given in Eq. (j2T))) . 
where now the M n in Eq. (TTJJ are given by the expres- 
sions in Eq. (|38[) . To define a point distribution in which 
the leading term in this expression vanishes, one can pro- 
ceed as follows. First one determines the location of the 
center of mass and the quadrupole moment of the elemen- 
tary tile. From the latter one can then find the principal 
axes, in which it is diagonal. In each tile of the tiling one 
then places on each such axis a pair of points, with their 
center of mass at that of the tile, at the appropriate dis- 
tance to produce the component of the second moment 
along the corresponding axis. The leading contribution 
to 5(k) should then be determined by the small-fc be- 
havior of 5 3 ' 3 (k), i.e., by the correlation properties of the 
third moment M 3 as in Eq. (|38|) . Because of the inversion 
symmetry in the point distribution in each tile, the first 
term in Eq. (|3"8|) vanishes. It is therefore the correlation 
properties of the third moment of the tiles alone which 
determines the coefficient of the term at order fc 6 . Given 
the results discussed above of the numerical studies 
for the quaquaversal tiling with a single point at the cen- 
ter of mass, we expect that such correlations are short- 
range, or at most very weakly long-range. One would 
expect therefore to obtain an S(k) cx fc 7 with 7 « 6. 

The generalization of this algorithm to higher orders 
is, in principle, straightforward (albeit evidently cumber- 
some as the order increases). Further one can seek to de- 
termine the minimal number of points per tile required 
to make the desired number of terms in Eq. (|20[) van- 
ish. Indeed the specific algorithm described above uses 
2d points, while it is easy to see that one needs only a 
smaller number to make the leading term in Eq. (|20[) van- 
ish: in d = 2, for example, it suffices to have three (rather 
than four) points to represent the quadrupole moment of 
the elementary triangular tile of the pinwheel tiling. 



IV. DISCUSSION AND CONCLUSIONS 

We have studied a class of point processes generated 
by placing one or more points in each tile of a tiling of 
M. d . We have assumed that the tiles have equal volume 
as this is expected to lead to the suppression of fluctua- 
tions at large scales characteristic of superhomogeneous 
point processes. We have shown explicitly that one can 
build superhomogeneous point processes with an arbi- 
trarily large exponent characterizing the small-fc behav- 
ior of the structure factor, and we have presented various 
examples. To our knowledge exact constructions of such 
point processes for the case 7 > 4 have not previously 
been given in the literature. 
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We have shown how in these algorithms the exponent 
7 depends (i) on the arrangements of the points ascribed 
to the tiles in the algorithm, and (ii) on the correlation 
properties of the shapes and orientations of the tiles. For 
the specific case of regular lattice tilings a non-trivial con- 
tribution to the small-fc behavior of 5(k) arises only from 
the former, with the coefficients in the small-fc expansion 
depending explicitly only on the variance of the (Senso- 
rial) mass moments of the points assigned to each cell. 
By arranging a sufficient number of points in a way which 
makes this variance zero for the first n moments, one can 
obtain 7 > 2n. In the case of an irregular tiling — we 
have considered the example of pinwhcel and quaquaver- 
sal tilings — an identical arrangement of the points in 
each tile can be sufficient to produce continuous SF and 
translational and rotational invariant superhomogeneous 
point processes. The exponent 7 then encodes informa- 
tion about the correlation properties of the shapes and 
orientations of the tiles. If these are short-range corre- 
lated, one obtains 7 = 4 placing a single point at the 
center of mass of each tile, and 7 = 2(n+ 1) if one places 
a number of points with all mass moments up to the n-th 
equal to that of the elementary tile. If there is, on the 
hand, long-range correlation in the shapes and orienta- 
tions of tiles, the exponent obtained will be modified in 
a way which depends on the nature of this correlation. 

Our results shed light on the meaning of the exponent 
7 > characterizing a superhomogeneous point process. 
Up to the value 7 = 4 previous explicit constructions of 
discrete processes (see, e.g., [ljl [26[) have shown that 
the increase of 7 can be associated with a suppression 
of fluctuations at large scales. Here we have seen that 
values 7 > 4 correspond indeed to an increased order in 
the arrangement of the points, but now at small scales: it 
is by changing how points are arranged within each tile, 
i.e., below a finite length scale (but subject always to the 
global constraints on fluctuations imposed by the tiling) , 
that we can increase the exponent. Thus to "undo" the 
order represented by an exponent 7 > 4 with respect to a 
system with 7 = 4 requires only the rearrangement of the 
system at small scales, while to "undo" that in a system 
with 7 < 4 require a coherent rearrangement of points 
on arbitrarily large scales (i.e. on scales inverse to the 
wavenumber range in which the exponent is measured). 

We have mentioned that our results are relevant in 
cosmology. Firstly the analysis given here makes more 
rigorous certain heuristic arguments used in this context 
regarding "causal constraints" on the generation of fluc- 
tuations from a uniform background [12l . [l3| . An algo- 
rithm like that described here, for the case of a single 
point placed at the center of mass of each tile, has been 
considered [l3| as a toy model for the generation of fluc- 
tuations starting from an exactly uniform mass density, 
by a physical process which conserves mass and momen- 
tum locally. The result 7 — 4 is obtained by assuming 
that space is divided into finite cells whose positions are 
uncorrelated. The latter assumption is in fact not consis- 
tent: the division of space into equal volume cells implies 



that their positions are necessarily correlated. Our more 
rigorous analysis shows that this exponent 7 = 4 does re- 
sult generically, however, if the shapes and orientations 
of these cells (i.e. tiles) are short-range correlated, i.e., 
have integrable correlation functions. 

Secondly, the generation of very uniform point pro- 
cesses is of relevance to the generation of initial condi- 
tions for numerical simulations of structure formation in 
the universe. In this context, to represent a given set 
of initial conditions, one must perturb appropriately (see 
[I?} for a detailed discussion) a point distribution rep- 
resenting as well as possible the uniform (unperturbed) 
universe. To understand the effects coming from this 
chosen point distribution (which are non-physical) it is 
desirable to have different algorithms which can generate 
such configurations. It is for this reason that a special 
case of the algorithm we have studied here has been built 
explicitly and studied numerically in this context [Io| . 
The analytical results we have given here complement 
these studies and give further algorithms for producing 
even more uniform point processes which may be useful 
in this context. We note again in this respect that ex- 
plicit algorithms for producing 7 > 4 have not previously 
been given. Such distributions in themselves provide in- 
teresting initial conditions (without any perturbation) for 
gravitational clustering, which have not previously been 
studied. 

We conclude with some further remarks on our results 
and some other directions for further work: 

• In our constructions of point processes we have al- 
ways constrained the center of mass of the points 
in each tile to coincide with that of the tile. We 
have done so because our goal here has been to 
generate point processes which are as uniform as 
possible. It is a simple exercise to redo our calcu- 
lation leading to Eqs. (| 1 5|) and (fTB)) when this con- 
straint is relaxed, i.e., allowing the center of mass 
of the particles (or particle) in each tile to be dis- 
placed randomly from that of the tile. The result 
is that the leading term in Eqs. (fT5|) is now at or- 
der fc 2 rather than fc 4 . For short-range correlated 
tilings the leading behavior of the SF will then be 
proportional to fc 2 . This result will be valid if the 
displacements of the center of mass of the particles 
within a tile with respect to the center of mass of 
the associated tiles have a finite variance. On the 
other hand, if the variance of these displacements 
diverges, the small-fc behavior of the SF is given 
by S(k) oc fc 7 where < 7 < 2, the value of the 
exponent depending on the precise behavior dis- 
placements PDF for large arguments. A detailed 
calculation of these cases for a randomly perturbed 
lattice can be found in Ref. [Ill, HB| ■ 

• While we have shown analytically the existence of 
point processes with arbitrarily large exponents 7, 
we have not done so for a case which is statistically 
translation and rotation invariant. In the latter 
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case our results for the exponent are expressed in 
terms of the small-fc behavior of the g n,m {k) which 
encode, as we have explained, information about 
the correlation properties of the shapes and ori- 
entations of the tiles. For the one such example 
which has been numerically studied (in Ref. [1 01 ] . 
the quaquaversal tiling with a single point at the 
center of mass) the result indicates an asymptotic 
behavior close to (but, as noted in Ref. (23|, slightly 
different to) that which would arise from a purely 
short-range correlation of the orientations of the 
tiles. Further numerical and analytical study of 
these points processes would clearly be of interest, 
in particular of the simpler pinwheel tiling. 

• We have made in our derivations here an assump- 
tion of analyticity at k = of the window function 
of the tiles, which corresponds to all moments of 
their mass distribution being finite. This is cer- 
tainly valid if the tiles are of finite extent. It may, 
however, include other cases which might be of in- 
terest, e.g., in d > 1 one may envisage that there is 
a non-trivial distribution of the shapes of the equal 
volume tiles, in which the extent of a tile is not lim- 
ited. One could also consider relaxing the assump- 
tion that the volumes of the tiles are strictly equal, 
admitting a distribution of volumes with specified 
correlation properties. In analogy with what has 
been found in certain algorithms for 7 < 4 [lfl Hi| , 
one would expect that such modifications would al- 
low the generation of point processes with leading 
non-analytic behavior, and any value of the expo- 
nent 7. Indeed one would expect non-analytic ex- 
ponents to be related either to the divergence of 
moments of such a distribution of extent or volume 
or to the presence of long-range tile-tile correla- 
tions. 

• While all our explicit examples have employed 
tilings which are congruent, our results for the 
small-A: behavior of the SF S(k) all apply only on 
the much weaker assumption of equal volume of 
the tiles. Thus for example we can apply these re- 
sults to any tiling generated by a deformation of 
tiles which leaves their volume fixed, which could 
encompass a large range of systems of physical in- 
terest (cells, foams, etc.). We recall in this respect, 
as remarked above, that a point placed randomly 
in each cell, rather than at the center of mass, leads 
to the restoration of the order k 2 term in the ex- 
pression given in Eqs. (fT5)) . 

• It is useful to briefly remark on the physical real- 
izability of superhomogencous point distributions 
with arbitrary positive but bounded values of 7. 
Such distributions are disordered to some degree 
and, although they are unusual, can be physically 
constructed. For example, the maximally random 
jammed (MRJ) state in three dimensions [28| is a 



special disordered sphere packing that can be re- 
garded to be a prototypical glass because it is per- 
fectly rigid and yet is maximally disordered. It is 
a superhomogeneous point distribution character- 
ized by an exponent 7 = 1 [6], but it is inherently 
a system out of equilibrium. While we have ex- 
amples like the equilibrium one-component plasma 
that has 7 = 2, can one devise equlibrium super- 
homogeneous point distributions in which 7 is ar- 
bitrarily large? The answer is apparently in the 
affirmative but it requires more than just pair in- 
teractions, namely, two-, three- and four-body in- 
teractions as shown in Ref. [l4j |. 
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APPENDIX A: FROM GAUSSIAN TO SPIN 
FIELDS 

In this appendix we show explicitly how to generate a 
regular lattice spin configuration with a given two-point 
correlation function. Such a configuration has been used 
as the starting point in the RBR algorithm described in 
Sect. HTB31 

The algorithm we propose is based on a mapping be- 
tween a set of correlated Gaussian variables {x(R)} with 
zero mean, and the spin set {s(R)} (where R is, as above, 
the generic lattice vector) . We do so because to generate 
a lattice set of correlated Gaussian variables with any 
possible desired correlation function is very simple. 

We denote by 



and 



c(R) = (a;(Ro)x(Ro + R)) 



X (R) = (s(R )s(Ro + R)) 



(Al) 



the two-point correlation functions of the Gaussian and 
the spin sets, respectively. We have used here the statisti- 
cal lattice translational invariance. Both c(R) and x(R) 
must have non-negative FTs as required by the Khint- 
chine theorem for stochastic processes (see e.g. Ref. [111). 

The starting point is the two variable joint PDF for 
correlated and monovariate Gaussian variables. Denot- 
ing by xi and xi two Gaussian variables at two lattice 
sites separated by the vector R, we have 



1 



2^ 4 -c2(R) 

a 2 (x\ + x\ — 2c(R)a;ia;2 



x exp 



2[ ( r 4 -c 2 (R) 



(A2) 
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where c(R) is defined in Eq. (|A1[) , and cr 2 = c(0) is the 
common variance of the Gaussian variables. 

The mapping we consider is the simplest possible: at 
the site R we fix s(R) = 1 if x(R) > and s(R) = -1 
otherwise. We want now to find the relation between 
x(R) and c(R). This can be done simply by noting that 
we can write 

s(R) = 26>[x(R)] - 1 , 

where 0{x) is the usual Heaviside step function. There- 
fore we can write 



x(R) 



dx\ dx2 



x exp 



to 2 v /<r 4 -c 2 (R) 
cr 2 (a; 2 + .Tj - 2c(R)xix 2 



1 .(A3) 



2[ct 4 - c 2 (R)] 

Performing the double change of integration variables 



2/i = xi 

yi=x 2 - [c(R)/a 2 ]a;i 



it is simple to rewrite Eq. (|A3|) as 



x(R) = 2 f 

Jo 



-I/?/(2<r a ), 



V2 



7T(J 



:f(A yi ) , (A4) 



where A = 



c(R) 



r^I^-C^R)] 



and 



erf(x) 



dt e 



is the usual error function. We now use the known equal- 
ity: 



dxe px erf(qx) 

which implies finally that 
2 



x(R) 



-arctg 



arctg ( -^z ) 



c(R) 



v^ 4 -c 2 (R) 



(A5) 



Therefore, given a lattice set of monovariate correlated 
Gaussian variables, we can map it onto a lattice set 
of spin variables with (s(R)) = and x(R) given by 
Eq. (£5). 



Asymptotics 



From Eq. (|A5j) it is simple to verify that 



X(0) = (a 2 ) = l. 



Moreover, as for R = |R| — > oo the correlation function 
c(R) must vanish, it is simple to verify that for suffi- 
ciently large R we have 



x(R) 



2 c(R) 



, i.e., x(R) an( i C (R) have the same asymptotic behavior. 
In particular if the Gaussian variables are long-/short- 
range correlated the spin variables are also long-/short- 
range correlated with the same scaling behavior. 

We can also give the condition of superhomogeneity for 
the spin lattice set. For the spin system this condition is 
simply 

R 

which gives the following more complicated relation for 
the correlation function of the Gaussian variables: 



x - 2 



R 



c(R) 



> 4 -c 2 (R) 



= 0. 
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